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The present status of the physics program, which led to the development 
of the Monte Carlo event generator PHOKHARA is described. The possi- 
bility of using the radiative return method in various aspects of hadronic 
physics, from the measurement of the hadronic cross section, to detailed 
investigations of the hadronic dynamics is emphasized. New results are pre- 
sented showing how to measure baryon form factors using the knowledge 
of their spin in baryon-antibaryon production with subsequent decay. 

PACS numbers: 13.40.Ks,13.66.Bc 

1. Introduction 

Tlie radiative return method proposed for tlie first time to serve as a 
tool in the hadronic cross section measurement [T] requires Monte Carlo 
event generator (s) as theoretical input . The method relies on the factor- 
ization properties of the differential cross section with photon(s) emitted 
from the initial states (ISR) and the possibility to solve problems caused 
by the photons emitted from the final hadrons (FSR). The first tools to 
meet the expectations from experimental groups were developed (EVA [2], 
EVA4pi [3]) using the structure functions method. That was a limitation, 
as not only the accuracy of the codes was not adequate to fulfill the ex- 
perimental expectations, but also within that framework more theoretical 
work is required to deal with possible double counting of the configurations, 
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where the hard photon is emitted at low angles. The adopted approach in 
the continuation of the research program, mainly fixed order exact calcula- 
tions, resulted in the development of the state-of-the-art Monte Carlo event 
generator PHOKHARAfT, 13 El El EJ [9], successfully used by the experi- 
mental groups BaBar, BELLE and KLOE working at meson factories. The 
hadronic cross section extracted using the radiative return method from 
the cross section of the reaction with emitted photons gives complemen- 
tary information to the scan method, traditionally used for such purposes. 
It is needed to calculate, via dispersion relations, the anomalous magnetic 
moment of the muon (for recent reviews see [10]) and the running of the 
electromagnetic coupling (for recent review see Profiting from the 

huge luminosities of the meson factories one can get valuable physical infor- 
mation, without building new accelerators and with accuracy competitive 
to traditional methods, as demonstrated [15] by the KLOE collaboration. 
As the method allows for the extraction of the hadronic cross section for 
energies from a production threshold to (almost) the nominal energy of the 
experiment, the B-factories have an access to data in the energy regions not 
covered previously by other experiments. Already now many new results 
exist, replacing the old measurements and/or covering new energy regions, 
crucial for a precise evaluation of {g — 2)^ and aqEoiQ'^)- The method orig- 
inally developed for the hadronic cross section measurement has, however, 
broader applications and can be used to study the hadron dynamics. Work 
started in this direction in [7], where it was shown how to extract nucleon 
form factors. It was continued in [8], where a method to test various models 
of the radiative (j) decays was proposed. Working along this line, a newly 
published analysis [13] shows how to use spin information on the decaying 
baryons to measure phase differences of their form factors. 



The Monte Carlo event generator PHOKHARA 6.0 was tested inten- 
sively at each stage of its development and the proved technical precision of 
the code is at the level of a small fraction (0.1-0.2) of a per mill. The pre- 
cision of the theoretical formulae used in the code is currently about 0.5% 
as far as the ISR is concerned. That accuracy is however not good enough 
to fully profit from more then 2 fb'^ data collected by KLOE [T5] at 
DAPHNE and further work is required to meet the growing expectations. 



The paper starts with a short description of the radiative return method 
in Section [2l The present status of the physics program for the precision 
hadronic physics with the PHOKHARA Monte Carlo generator is outlined 
in Section [3l A short summary is presented in Section [H 
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2. The radiative return method in short 

The radiative return method reUes on factorization properties of the 
cross section with photon emitted from the initial electron or positron 

da{e~^e~ hadrons + 7(ISR)) = 

H{Q^,e^) do-(e+e- ^ hadrons, s = Q^) , (1) 

where is the invariant mass of the hadrons. A similar factorization holds 
also, when more than one photon is emitted from initial states. The function 
H{Q'^,9j) (or a more complicated function for multi-photon emission), at 
relatively low energies of meson factories, is given with high accuracy by 
QED only and is thus well known. From this follows that a measurement 
of the differential (in Q"^) cross section of the reaction e^e" hadrons + 
photons allows for a cross section a{e^e~ hadrons) extraction for energies 
from the production threshold to almost the nominal energy of a given 
experiment. 

The presence of the contributions from photon(s) emitted from final 
hadrons has to be treated as a background. It has to be studied carefully as 
the models of photon emission from hadrons are not well established. The 
region of the hadron invariant masses, which is of main interest, is below 3 
GeV and thus the role of the FSR at the <j) factory DAPHNE is far more 
important then at B- factories (BELLE and BaBar). The reason is purely 
kinematical. To obtain a low invariant mass of the hadronic system one 
needs a very energetic photon emission at B- factories, which is not the 
case for the relatively low energy of (p- factory. The typical kinematical 
configuration at a B- factory is thus an energetic photon and hadrons going 
in opposite directions. As the FSR contributions are enhanced only in the 
kinematical regions, where the hadrons and photons directions overlap, the 
FSR is naturally suppressed at B-factories and it is not suppressed at the 
(p- factory DAPHNE, where special care is needed to deal with it. This is 
shown in Fig. [H where the relative contribution of the leading order FSR 
corrections is plotted for the e^e~ — > Tr~^TT~j differential cross section. Even 
if the FSR contribution at DAPHNE is sizable, when no event selection is 
used, it is relatively easy to choose an event selection, which suppresses 
this unwanted background. A simple choice is also shown in Fig. [H but 
more sophisticated (and efficient) solutions can be found, as it was done 
by the KLOE collaboration in their pion form factor measurement [12]. As 
the leading order contributions are suppressed at B- factories by the choice 
of specific kinematical configurations, it is not a surprise that the next to 
leading corrections, which are not suppressed by that choice, are bigger then 
the leading corrections. They amount up to a few percent and do depend 
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Fig. 1. The relative leading order FSR contribution to the differential, in the in- 
variant mass of the pion pair (Q^), cross section of the reaction e'^e~ tt^tt^^ 
for DAPHNE and B- factories energies. 



on the invariant mass of the pion pair, and therefore have to be taken into 
account, when aiming for a precise measurement. 

For DAPHNE energies the event selection used by KLOE to suppress 
the FSR contributions does not ahow for a measurement of the pion form 
factor in the threshold region. That region might be important, as the 
contributions from low invariant masses are enhanced by the kernel func- 
tion in the dispersion relations for the muon anomalous magnetic moment. 
Releasing the cuts allows for measuring in the threshold region, but in the 
same time the model dependence of the FSR contribution becomes a serious 
problem [IB]. Detailed studies of the FSR, also for the specific case of a (j)- 
factory, where the radiative cj)- decays play a role, were performed in many 
theoretical papers El El [T7] and that discussion will not be repeated here. 
One comment is, however, in order: the final solution to that problem can 
only come through close collaboration between theoretical and experimental 
groups and the analysis of both on peak and off peak data is crucial for the 
success of the program [16] . 

3. An overview of the present status of the theory research 
program and new challenges 

The current version of the PHOKHARA 6.0 Monte Carlo event genera- 
tor is a product of many years of theoretical investigations, calculations and 
code testing. It relies on virtual (-|-soft) radiative corrections to the ISR 
calculated in [18] and ISR hard photon corrections calculated by means of 
the helicity amplitudes and implemented into an efficient Monte Carlo event 
generator in [31[5]. The FSR at the next to leading order was investigated, 
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and implemented into the event generator, for muon and pion pair produc- 
tion in [6l[9]. In addition, very specific contributions to FSR, important at 
the (f)- factory DAPHNE and coming from radiative (p- decays, was stud- 
ied in [8]. All the implemented parts of the developed code are tested to 
achieve a relative technical precision of few times 10~^. Independent tests 
were performed in [19], where the authors compared ISR contributions of 
muons present in PHOKHARA against the KKMC event generator [20] . An 
excellent agreement between non exponentiated version of the KKMC and 
PHOKHARA was found, while higher order corrections, not implemented 
yet in the PHOKHARA, give at most two per mill contributions for the 
invariant masses of the muons relevant for the radiative return method. 

With PHOKHARA 6.0 one can now generate the following final states: 
Tr+vr", K+K-, K'^K^, pp, nn, 7r+7r"7r°, 27r+27r-, 27r°7r+7r-, A(^ 

7r+p)A(^ TT^p), accompanied by one or two ISR photons. The FSR cor- 
rections are implemented only for 7r"'"7r~, jj.'^ and K^K~ ^ while in the 
prepared new release they are implemented also for pp and 7r"'"7r~7r^ final 
states. The narrow resonance {J/ip and ^/)(2S')) contributions to 7r"'"7r~, 
fi~^fi~, K~^K~, K^K^ will also be implemented there. The A pair produc- 
tion and decays are implemented at the leading order only [13j, but as the 
expected number of events is modest, the accuracy of the code is sufficient 
for the description of this process. The spin asymmetries and spin-spin 
correlations of the lambdas provide information about real and imaginary 
parts of the lambda form factors. To measure them, one needs only the 
information on the angular distributions of the produced lambdas and pi- 
ous coming from their decays [13j . A step towards such a measurement was 
done by BaBar in [21j, but only a limited part of the information contained 
in the data was actually used. 

The model describing the four pion channels was improved, based on 
experimental information from BaBar [22] and CMD2 |23tl24j charged mode 
measurements together with CLEO [25j and ALEPH [26] spectral functions. 
The preliminary results indicate some isospin symmetry breaking effects, if 
one also adds information from the preliminary BaBar measurement of the 
neutral mode [27] . 

Even if relatively few final states are implemented in the distributed version 
of PHOKHARA, as compared to the plethora of available final states, the 
implementation of the missing channels by a potential user is not difficult, 
at least for the ISR part, due to the modular structure of the program. This 
was done for example by the BELLE collaboration and used in [28]. 

The 0.5% accuracy of the ISR corrections in PHOKHARA will soon 
be the biggest contribution to the error in the pion form factor extraction 
by KLOE. Inclusion of the leading logarithmic corrections from the second 
loop to one photon emission, leading logarithmic corrections at one loop 
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level to two photon emission and implementation of the three hard photon 
emission is expected to bring the accuracy to the 0.2% level and this will 
be the priority of the physics program of the group in the near future. It 
will follow the accuracy improvement [29] of the BABAYAGA code used by 
KLOE for the luminosity measurement. 



4. Summary 

The present status of the radiative return research program was outlined 
and plans for the near future work towards further improvements of the 
PHOKHARA Monte Carlo generator were sketched. 

The publication is based in a big part on results obtained in collabora- 
tion with J. H. Kiihn, E. Nowak-Kubat and G. Rodrigo. The authors are 
grateful for many useful discussions concerning experimental aspects of the 
radiative return method to members of the KLOE, BaBar and BELLE col- 
laborations, mainly Cesare Bini, Achim Denig, Simon Eidelman, Wolfgang 
Kluge, Debora Leone, Stefan Miiller, Federico Nguyen, Evgeni Solodov, 
Graziano Venanzoni and Ping Wang. 
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